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Abstract 

Optimal two-dimensional (2D), three-dimensional (3D) wave launching configurations are pro- 
posed for enhanced acceleration of charged particles in magnetized plasmas. A primary wave 
is launched obliquely with respect to the magnetic field and a secondary, low amplitude, wave is 
launched perpendicularly. The effect of both the launching angle of the primary wave, and the pres- 
ence of the secondary wave is investigated. Theoretical predictions of the highest performances 
of the three-dimensional (3D) configurations are proposed using a Resonance Moments Method 
(RMM) based on estimates for the moments of the velocity distribution function calculated inside 
the resonance layers (RL). They suggest the existence of an optimal angle corresponding to non 
parallel launching. Direct statistical simulations show that it is possible to rise the mean electron 
velocity up to the order of magnitude as compared to the primary wave launching alone. It is a quite 
promising result because the amplitude of the secondary wave is ten times lower the one of the first 
wave. The parameters used are related to magnetic plasma fusion experiments in electron cyclotron 
resonance heating and electron acceleration in planetary ionospheres and magnetospheres. 

PACS numbers: 52.35.Mw, 94.20.Rr, 94.30.Hn, 52.25.Xz 
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Electron acceleration due to external radio frequency waves in a strong magnetic field has 
long been recognized as an important effect in a wide variety of problems ranging from plasma 



heating and current drive in fusion devices 
radiation belts during geomagnetic storms 



probing 




to electron acceleration in the Earth's 
active ionospheric and magnetospheric 



It is well known that wave-particle interactions are most efficient when the particles are 
in resonance with the waves. The resonance conditions 
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define some regions 'of sensitivity' in wave-particle parameter space which can be described 
as Resonant Layers (RL). Here, N is the harmonic number, u = 27r/ is the wave frequency, 
k\\ and v\\ are the components of the wave vector and the electron velocity parallel to the 
constant magnetic field Bq, 7 = (1 — f^/c^)~^/^ is the relativistic factor and = eBo/nie the 
gyro-frequency. This view allows to develop significantly previous simulations of c harg ed par- 
ticle fluxes and plasma disturbances in ambient magnetic field in 12, 13, 14, 15| . There 
is a special case when electrons permanently staying in the RL (^. Such a phenomenon has 



been referred to as autoresonance l2G 
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22 and its conditions are known as the cyclotron 

auto-resonance maser (CARM) conditions^- 

Several mechanisms have been explored 

for maintaining the synchronization between electrons and waves not fully satisfying these 
CARM conditions such as changing the profile of the guide magnetic field or varying the 



wave phase velocity 
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25| . Recently, the use of two parallel counter-propagating waves has 



been considered l26l l27l . Numerical tests 28] have shown that the two-wave scheme may 
lead to higher averaged parallel velocity. The stochastic acceleration mechanism for elec- 
trons in a plane monochromatic electromagnetic wave propagating obliquely to the external 
magnetic field has also been studied It was found that it is easier to accelerate 

electrons to high energies with increasing propagation angle when the electron motion be- 
comes stochastic and the parallel phase velocity of wave is supraluminous (nii = knc/uj < 1). 



Furthermore, Karimabadi and Angelopoulos 



33| studied the interaction of charged particles 



with a packet of obliquely propagating plane monochromatic electromagnetic waves under 
the special condition (of equal cos 9i for all the waves, where for the z-wave - refraction 
index, 6i - the propagation angle). This condition allowed the system to be reduced to 
two degrees of freedom and the particles can be accelerated through a process of Arnol'd 
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diffusion 3J] in two dimensions. 

The majority of the existing works are based on the description of a single particle 
dynamics in one (or more but under condition of equal cos6'i) plane monochromatic radio 
frequency waves. 

The coherent acceleration of nonrelativistic ions interacting with two and more elec- 
t..o.tat,c waves a ..fo™ ™et,c field has been studied by -ent.y by Q. as a 
generalization of analysis of 0, Q] by including wavenumbers along the external magnetic 
field. 

In this paper, a mechanism is discussed for the acceleration of electron populations re- 
sulting from the effect of crossing electromagnetic waves propagating in a dispersive medium 
according to the geometry represented in Figure Q (the condition of equal rii cos 6i for the 
two waves is thus clearly broken). To analyze this mechanism, the resonance moments (RM) 
of the distribution, i.e. velocity moments computed in the RL only, are evaluated. The first 
order RM suggests that a peculiar 6 results in a maximal averaged parallel flux. Although 
the RM approach has to be considered as an approximation, this prediction is reasonably 
confirmed by direct statistical simulations. Moreover, the two-wave scheme allows to rise 
the mean electron velocity up to one order of magnitude when compared to the one-wave 
scheme, based on the primary wave only. 



X 
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FIG. 1: Schematic picture of the electro- magnetic configuration. 



The electromagnetic configuration that is considered (Figure is the combination of a 
strong magnetic field (assumed to be along the z direction), a primary wave propagating 
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obliquely with respect to the magnetic field and a secondary wave propagating perpendicu- 
larly to the magnetic field. As a first step, to simplify direct particle simulations both the 
primary and the secondary waves are assumed to be in the plane {x, z). This simplification 
can give not so impressive effect as compared to 3D launching. Nevertheless, it will give 
a first estimate of the secondary wave effect and motivation to develop more realistic 3D 
launching code that will be closer to real experimental setups. 

This electromagnetic configuration is not an attempt to satisfy the resonance condi- 
tion during a long time being close to the cyclotron auto-resonance maser (CARM) condi- 
tions [3]. Rather, a large population of electrons is considered and only the average effect 
of the waves on the population is considered. The fraction of electrons that are close to 
the condition (Q), that corresponds to a resonance layer (RL) in the velocity space, becomes 
then as important as the time these electrons remain resonant. The secondary wave does not 
carry any parallel momentum and cannot induce any net parallel motion of the electrons. 
The purpose of this secondary wave is to maintain a pseudo-equilibrium velocity distribution 
in which the RL is continuously re-filled. Indeed, the combined effect of the two-waves and 
the magnetic field yields a stochastic motion during which the synchronization between the 
waves and the gyro-motion of the electron is, on average, more favorable for transferring 
momentum to the electrons. 

The description of the electron trajectory (r, p) requires a relativistic treatment and is 
derived from the time dependent Hamiltonian: 



and the trajectories of the electrons are determined by the initial conditions and by the 
Hamilton equations. This picture corresponds to a test-electron in an external electromag- 
netic field. Assuming that the electrons interact with two monochromatic waves propagating 
in cold plasma at angles 6i and 62 with respect to the guide magnetic field, the total vector 
potential can be written as follows: 



where ei and 62 are two unit vectors in the plane {x — z). These vectors, as well as the an- 
gles ipi and ^/'2 and the refraction indices ni = ||ki||c/a;i and n2 = ||k2||c/a;2, are determined 
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H = ^/rri^c^ c2(p -h eA)^ , 



(2) 



A 



+ — 




(3) 



by the Appleton-Hartree dispersion relation in the cold homogeneous plasma approxima- 
tion 35] . Collisions with other particles, electrons or ions, are neglected. The dynamics of 
the electrons in the electromagnetic configuration Q is chaotic and unpredictable analyti- 
cally in general. Electrons with slightly different initial position and velocity may experience 
drastically different evolutions. The exact analytical prediction of the average effect of the 
waves is thus out of reach. It is however possible to anticipate the existence of a RL (HJ 
in the velocity space to estimate the possible net effect of the waves by computing the RM 
defined by: 



for any function g of the velocity. It represents the density of electrons close to the RL for g = 
1 and the mean parallel fiux of the electrons on the layer for g = v\\. These quantities should 
give some estimate on the efficiency of the electron-wave interaction. The RL corresponds to 
an ellipse in the velocity space and the integral / can be evaluated analytically using elliptical 
coordinates. For instance, assuming a Maxwellian distribution, /(v) = Cexp(— /?v^), the 
evaluation of / for (7 = fy, which will be denoted /* can be done explicitly. This quantity 
corresponds to the averaged parallel fiux of the particles that belongs to the RL. It is, 
at least indirectly, related to the net averaged parallel current produced by the electron- 
wave interactions. Indeed, these interactions tend to remove electrons from the RL while 
the pseudo thermalization of the electrons due to the combined effect of the two waves is to 
refill constantly the layer. The thermalization is thus expected to add a net averaged parallel 
velocity proportional to /*. Of course, /* gives only a rough indication of the efficiency of 
electron-wave interactions and the final averaged velocity cannot be deduced directly from 
it. The exact expression for I* is quite long. It is thus more illustrative to present the 
Figure 121 in which /* is a function of 61, the angle of the primary wave. The dimensionless 
wave amplitudes and the quadratic plasma frequency are defined by = Ai^2^/cBq, 
Co = (cUpe/fi)^- The other parameters correspond to the values used in the simulations 
described below. 

As expected, for perpendicular propagation 9i = 90°, no averaged parallel fiux is observed. 
For two sets of parameters, there is a clear maximum of the averaged parallel fiux of the 
particles in the RL for 61 ^ 0°. The explanation for such a phenomenon is that the averaged 
parallel velocity induced by the electron-wave interaction depends on both the angle of 
propagation of the primary wave and the number of electrons that are close to the RL. 




v£ res. layer 




(4) 
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FIG. 2: Averaged parallel flux of particles in the RL in arbitrary units for a Maxwellian distribution 
for increasing densities, cq = 0.3 solid line, cq = 0.6 dashed line and cq = 1.99 dotted line. 
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FIG. 3: Averaged electron velocity versus 6i (in degrees) for the one- wave scheme with = 0.1, 
^2 = (dotted lines) and the two-wave-scheme with Ai = 0.1, A2 = 0.01 (solid lines) at time 
n t = 7000, for low-density runs eo = 0.3. 

Smaller angles of propagation correspond to higher parallel momentum carried on by the 
wave. However, the RL condition is compatible with larger numbers of electrons for higher 
angles of propagation (at least assuming a Maxwellian distribution). It should be noted that 
here /* has been computed using only the fundamental N = 1 RL. Contributions from the 



6 



higher harmonic layers decrease rapidly because these layers are more and more symmetric 
and because the absolute value of the resonant velocity increases towards high values for 
which the electron density is very small. 
70 
60 
50 
40 
30 
20 
10 

-0.2 0.0 0.2 0.4 0.6 0.8 V^/ C 

70 

60 
50 
40 
30 
20 
10 

/ 
-0.2 0.0 0.2 0.4 0.6 0.8 Mj C 

FIG. 4: Final probability distribution of the parallel velocity for 9 = 15° and Ai = 0.1 for the 
one- wave scheme with ^2 = (top) and the two- wave-scheme with A2 = 0.01 (bottom). The solid 
lines (eo = 0.3), dashed lines (eo = 0.6) and dotted lines (eo = 1.99) correspond to increasing 
densities. 

The Hamiltonian equations for the r and p are solved using a 4th order Runge-Kutta 
method. The time step is adapted to ensure that the solution of a redundant evolution 
equation for H remains close to the expression 0. The initial velocity distribution has a 
temperature of the order of 1 keV. A population of 5 IC^ electrons is used in each simulation. 
Although running the code with larger populations is not an issue, no further information 
is derived from these larger runs, except of course more converged statistics. 
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Numerical results confirm that the angle dependency is not trivial and that parallel prop- 
agation (^1 = 0°) for the primary wave is not always optimal ^29]. Three sets of simulations 
are presented hereafter. The parameters for these simulations are relevant in today toka- 
mak plasma. In particular, the primary wave corresponds to the second harmonic of the 
cyclotron frequency and the secondary wave to the third harmonic Right Hand Polarized 
modes which are frequently used for instance in the TCV tokamak experiments js^. The 
value of the constant magnetic field is 1.42 T for all simulations. Three electron densities 
have been considered: ~ 0.6 10^^ (cq = 0.3), rig ~ 1-2 10^^ (cq = 0.6) and 
He ~ 3.9 10^^ (eo = 1.99). The wave amplitudes are Ai = 0.1 and A2 = 0.01 in the 
two-wave scheme and ^2 = in the one-wave scheme. They correspond to power fluxes 
which are by orders of magnitude higher than that achievable on gyrotron used in today 
tokamak. However, preliminary computations using the Resonant Moment Method suggest 
three-dimensional electromagnetic wave configurations are very promising for larger accel- 
eration of charged particles in an external magnetic field with even lower wave amplitudes. 
In such a case the wave vectors and the magnetic field are not supposed to be co-planar and 
create a fully three dimensional system. On the other side, the required powers might be 
achievable by free electron maser (3.7| even for experiments with 2D launch configuration as 
predicted by our direct particle simulations. 

Also, the parameters Cq = {urpe/^)"^ = 0.1 — 0.3, correspond to the nighttime ionosphere 
at approximately 130km, A=0.1 - to a power flux 5W/cm'^, and a frequency 2.6MHz of the 
primary wave. These parameters are close to ones considered in [2^ 0] for single wave 
acceleration. 

Figures El show a significant increase of the average parallel velocity for cq = 0.3 due to 
the secondary wave. Moreover, the angle dependency of the average parallel velocity appears 
to be maximal in the range 61 = 10° — 60°. This range thus appears to significantly differ 
from the privileged value of the RM I* (Fig. 2). This is not too surprising since the RM 
have been computed assuming a Maxwellian distribution with zero mean. Thus, although 
the global picture from the RM description is reasonable, a more precise iterative approach, 
taking into account the averaged velocity suggested by the RM would be required for more 
accurate predictions. The probability distributions of parallel velocity (Figures E)) observed 
at the end of the simulations Q t = 2 10^ indicate that, in the two- wave scheme a much 
larger number of electrons have had the occasion to interact with the primary wave and 
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the distributions of velocity exhibit two well marked maxima. For cq = 1.99, the density is 
very close to the cut-off value of the wave propagation and almost no effect is observed in 
both the one-wave and the two-wave scheme. Also, if Ai is too small, no averaged parallel 
velocity is observed at all. Preliminary tests seem to reproduce the threshold previously 
observed 

0,0- 

This paper presents a mechanism for enhancing acceleration of a population of electrons 
using crossing electromagnetic waves propagating at different angles with respect to an 
external magnetic field in a dispersive medium. The existence of optimal angles of prop- 
agation for the primary wave is suggested using the evaluation of resonant moments and 
is confirmed by direct numerical simulations of the electron trajectories. A secondary low 
amplitude perpendicular wave is used to induce a stochastization of the electron trajectories 
and, consequently, to maintain a pseudo-equilibrium. Although measures of the distribu- 
tions (Figures H} clearly show a departure from a thermal equilibrium, the stochastization 
effect of the secondary wave yields a clear increase of the average parallel electron velocity. 
It is a quite promising result since the amplitude of the secondary wave is ten times lower 
the one of the first wave. 

The authors are grateful to Professor R. Balescu, Dr. B. Weyssow, Dr. I Pavlenko, and 
Dr. R. Kamendje for useful discussions. D.C is researcher of the Fonds National pour la 
Recherche Scientifique (Belgium). This work has been supported by the contract of associa- 
tion EURATOM - Belgian state. The content of the publication is the sole responsibility of 
the authors and it does not necessarily represent the views of the Commission or its services. 
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